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ABSTRACT 

Orthotropic plate theory has been increasingly used to 
model plate-stiffener combinations typical of those which 
are used in ship hulls. The behaviour of a thin plate with 
large deflections is described by two nonlinear partial 
differential equations of equilibrium and compatibility. 

The orthotropic form of these equations is derived and 
solved for a large aspect ratio plate with long edges simply 
supported and short edges fixed as boundary conditions. 
Buckling and post-buckling regions are investigated under 
combined loadings of lateral pressure, inplane edge com- 
pression and edge shear. Results are presented for virtual 
aspect ratios 1/1.5 to 1/6 and both isotropic and ortho- 
tropic plate properties in the form of design and behaviour 
charts . 
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' NOMENCLATURE 


y= 


ordinate in long direction 


X 


= ordinate in short direction 


a 


= plate length in x-direction 


b 


= plate breadth in y-direction 


3 


= — = aspect ratio 

3. 


U, V, w 


= displacements of a point in x-, y-, and z- 
directions respectively 


e , e , e 
x y xy 


= middle-plane strains 


D , D 
x y 


= flexural rigidity of orthotropic plate in x- 
or y-directions respectively 


D xy 


= effective torsional rigidity of orthotropic 
plate 


V , V 
x y 


= Poisson's ratio of orthotropic plate in x- or 
y-directions respectively 


E , E 
x y 


= modulus of elasticity in x- or y-directions 
respectively 


G 


= modulus of elasticity in shear 


N , N , N 
x y xy 


= middle-plane loads per unit length 


M , M 
x y 


= bending moment in orthotropic plate acting 
around a line perpendicular to x- or y-axis 
respectively, per unit width 


M 

xy 


= twisting moment in orthotropic plate 


b 


= nondimens ional deflection coefficient 



mn 
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F = 


Airy's stress function 


J = 

X 


1 

hE 

Y 


h = 


plate thickness, orthotropic 


J = 
Y 


1 

hE 

X 


2J 

xy 


1 T T 

Gh x y y x 


2d 

xy 


D v + D v + 4C 
x y y x 


C = 


Gh 3 

12 


e = 


■ 12 ^r (1 ~ v x v y ) 


p , p = 

x y 


= total loads in x- and y-directions respectively 


q = 


= uniform lateral load, per unit area of plate 


s = 


= constant inplane shear load, per unit length 


p 8 


D ' 
cl V 

= £ = virtual aspect ratio of orthotropic plate 

v] X 


n = 


D 

= — t = torsion coefficient of orthotropic plate 

Id d 
si x y 


Y = 


J 

xy 

|j J 
\l x y 


N* = 
X 


N x a2 

= = nondimensional in-plane load x-direction 

h 2 D 

X 



12 



N* = 

y 


N b 2 

= — ^ — = nondimensional in-plane load 

TT 2D y 

y-direction 


Q* = 


Qb ** 

= — 2 = nondimensional lateral load 

^hD y 


s* = 


_ __^_ab_ _ ncn( ^i raens ^ ona ^ edge shear load 

TT 2 D 

xy 


b , a = 
e e 


= effective width, x- or y-directions 
respectively 


°x' °y ; °x' °y = 


= membrane stresses in x- and y- 

directions respectively; nondimensional 


T xy ' T xy 


= shear stress in xy plane; nondimen- 
sional 


a bx' °by ' °bx' °by 


= bending stresses in x- and y- 

directions respectively; nondimensional 


“T • rr "fa ~ 

T bxy ' T bxy 


= twisting stress in xy plane; non- 
dimensional 


°tx' a ty' T txy 


= total stresses 


a tx' a ty ' T txy 


= nondimensional total stresses 


a tx a ty' T txt 


= nondimensional total stresses 


°1 , 2 ' a l,2 


= maximum and minimum principal stresses, 
plane stress; nondimensional 


a x ' °y = 
e 2 e 


= edge membrane stress in the x- or y- 
direction respectively 



nondimensional bending moments 
per unit length 

nondimensional bending moments per 

unit length due to curvature in x- 

direction only (M* ) and y-direction 

x o 

only (M* ) 
coefficients 

Von Mises combined stress; non- 
dimensional 

moments of inertia of the stiffeners 
with effective plating in the x- and 
y-directions respectively 
moments of inertia of the effective 
plating alone in the x- and y- 
directions respectively 
spacings of the stiffeners extending 
in the x- and y-directions respec- 
tively 

equivalent thickness of the plate 
and the stiffeners (diffused) in the 
x- and y-directions respectively 
thickness of plate alone 
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1. Introduction 



Of major concern to the naval engineer is the reaction 
of ship bottom plates to a combination of inservice loadings. 
In general, the plate-stiffener combinations will be sub- 
jected to inplane compression or tension edge loads due to 
ship hogging or sagging, lateral hydrostatic loads and in- 
plane edge shear loads. With increasing compression and 
shear loads, a critical state will be reached where buckling 
will occur. Plate deflection will exceed plate thickness, 
however, the plate may still carry considerable loads. High 
slenderness ratio (b/h) plates are encountered in modern 
longitudinally framed ship designs. The yielding load may 
be considerably higher than the buckling load and investi- 
gation of the plate-stiffener residual strength in the post- 
buckling region is of great importance. 

Orthotropic plate theory has been increasingly applied 
to ship structures. With this analysis, the plate-stiffener 
combination is modeled as an equivalent flat plate with 
elastic properties that are different in two perpendicular 
directions. Schade [12, 13] and Mansour [9, 10]*, among 
others, have worked with this concept. 

The small deflection theory is useful only in the pre- 
buckling region where deflections of a plate are small in 



*Numbers in brackets designate references in Chapter 8. 
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comparison with its thickness, the assumption of no deforma- 
tion in the middle plane of the plate is made in this case. 

If deflections are not small, ^ > j, the strain in the middle 
plane is no longer negligible and must be considered. Von 
Kantian 1 s equations describe the behavior of an isotropic 
plate under combined lateral and inplane loads. The non- 
linear equations describe the behavior of the plate in both 
small and large deflection regions. 

Considerable application has been made of von Kantian' s 
equations. Levy [4, 6] investigated square and rectangular 
plates simply supported and subjected to both inplane edge 
loading and normal pressure. Coan [3] included the effects 
of small initial curvature for this case. Levy [5] , using 
two fixed and two simply supported ends, solved the equations 
with edge loading only for an aspect ratio of 4, and studied 
pure shear for the simply supported square plate [7, 8], 

The square plate with variations of the boundary conditions, 
subjected to edge loading only, was studied by Yamaki [17] . 
Payer [11] applied von Kantian’ s equations to deep web frames 
and included uniform edge shear in addition to compressive 
edge stress and normal pressure. Shultz [15] investigated 
wide plates having aspect ratios of $ = 1.5 to 8 for con- 
ditions of a transversely framed ship. All of these 
analyses apply only to isotropic plates. Mansour [9, 10] 
extended von Kantian' s equations to consider slightly rec- 
tangular orthotropic plates under various boundary, loading 
and initial deflection conditions. 
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A large aspect ratio orthotropic plate with the short 
edges fixed and the long edges simply supported is considered 
in this paper. The plate-stiffener combination is subjected 
to loads of inplane edge compression, uniform edge shear and 
uniform pressure normal to the plate. Figure 1 illustrates 
the loadings and coordinate system used. 

The orthotropic form of von Harman's equations are 
solved similar to the theoretical analysis of Mansour [9, 10]. 
The IBM System 370 model 165 is employed to produce the 
numerical results. "Design" charts of deflection, effective 
width and bending moment are given for orthotropic plate 
virtual aspect ratios of 1/1.5 to 1/6 with inplane edge 
compressive and lateral loadings. In addition, "behaviour" 
plots of plate centerline deflection and total principal 
stresses are given for virtual aspect ratios of 1/1.5 and 
1/1.25 for various combinations of inplane edge compression, 
edge shear and lateral loadings . 
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2. Formulation of the Problem 



2 . 1 Orthotropic Material Properties 



An orthotropic material has different elastic proper- 
ties in two perpendicular directions. For plane stress in 
the xy plane, the stress-strain relations are 



c 



x 



E 

l-v x v y (e x + v y e y ) 



E 

a = r: — — (e + v e ) 

Y 1-v x v y ^ x x 



( 1 ) 



T = Y G 
xy 'xy 



from energy symetry 



E v = E v 
x y y x 



which gives 4 independent elastic constants. 



2 . 2 The Rectangular Orthotropic Plate with Large Deflection 
(Figure 1) 

, , w 1 . 

For thin plates with large deflection, ^ a satis- 

factory approximate theory makes the following assumptions: 
1. Points initially on the normal to the middle plane 
of the plate remain on the normal after bending. 
(Disregard of shear deformation.) 
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2. The normal stresses in the direction transverse to 
the plate can be disregarded. (Plane stress where 



3. Hooke's Law relating stress-strain applies. 



For the large deflection theory, deformation in the 
middle plane of the plate must be considered. These strain 
components therefore include the effect of deflection and 
are approximately 



e 

x 



3u 2 

3x 2 K dx 



£ 



y 



3v ^.SW, 2 
3y 2 K dy } 



(2) 



Y 



xy 



3u 3v_ + 3w . 3w 

3y 3x 3x * 3y 



Equilibrium of forces and moments on an element of a 
plate produces the equilibrium equations 



and 



3N 



x 



3N 



3x 



= 



3y 



~3y "lx" 



= 0 



(3) 



3 2 M 



x 



3 2 M. 



3x" 



- 2 - 



xy; 



3x3y 



3 2 M 



dy‘ 



= -(q + N 



3 2 w 
X 3x 2 



+ 2N. 



xy3x3y 



9 * w + N liSL) 

y 3y 2 



(4) 
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Introducing the familiar Airy's stress function F which 
satisfies equations (3) 



N = 
x 9y 2 



N = ^ ; N 

y ax 2 



xy 



3 2 F 

9x9y 



(5) 



From substitution of large deflection strains (2) into 
the stress-strain relations (1) and using the definition of 
bending and twisting moments, the moment-curvature relations 
for an orthotropic material are derived 



M = -D + v 

x x y 



9j_W) 
9y 2 



M = -D + v 

y y 3 Y 2 x 



9 2 W j 

9x 2 




2C 



9 2 w 
9x9y 



where the rigidity coefficients are defined as 

E h 3 

D = x 

x 12 ( 1-v X v y) 

E h 3 

D = 1 

y " I2(i-v x v y ) 

and 



( 6 ) 



C 



_ Gh 3 
" 12 



Substitution of equations (6) and equations (5) into equation 
(4) , the "equilibrium" equation is obtained as 
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D + 2D 

X 3x“ Xy 



3 4 w 
3x 2 3y 5 



+ D 



y 



3 4 w 

3y 4 



= q + Li I . i!w 

3y 2 3x 2 



2-JiZ 

3x3y 



( 7 ) 

3 2 w 3 2 F 3 2 w 

3x3 * 3x* ' 3y* 



where 




D v + D v 
x y y x 



+ 4 C 



To obtain the "compatibility" equation, the strain equations 
(2) are first differentiated and combined to eliminate dis- 
placements u and v 



3 2 e 



x 



3 r 



3y‘ 



xy 



3x3y 



3 2 e, 



3x‘ 



= ( $2” ) 2 - 

dxdy 



3 2 w 
3x 2 



3 2 w 

3y 2 



( 8 ) 



Differentiation of equations (1) , substitution into equation 
(8) along with (5) and using the equilibrium equations (3) , 
the "compatibility" equation is obtained as 



where 



and 



J + 2J . ■ - 3 - 4 .F , + j ^ 

x 3x 4 X ^3x 2 3y 2 ^3y 4 



= ( 



3 2 w 
3x3y 



3 2 w 
3x 2 



3 2 w 

3y 2 




(9) 



- v J 
x y 



v 



J 

y x 



2J xy = tfe 
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Equations (7) and (9) are fourth-order non-linear 
partial-differential equations that describe both small and 
large deflections of an orthotropic plate. This includes the 
buckling and post-buckling behavior of the stiffened plate. 
Substitution of the isotropic material properties 

£ 

E = E = E;v = v =v; G = 0 . ■ . 
x y x y 2(l+v) 

result in the familiar von Kantian' s equations. 

Solution of these equations with 16 boundary conditions 
results in the two functions F and w. The boundary con- 
ditions are eight support conditions and eight edge load/ 
displacement conditions. 



2 . 3 Bending and Membrane Stress 

The membrane stresses are determined from equations (5) 
as N 



0 - TT 

x h 



1 3J\F 
h 3y 2 



N 



„ _y 1 9 f 

a = -rf- = s— 



(10) 



9x" 



_ N xy _ _ 1 3 2 F 
T xy h h 9x9y 



Noting that the maximum normal stress acts on those sections 
parallel to the xz or yz planes, and using equations (1), 

(2) and (6) the bending and shear stresses are obtained as 
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E M 

, x x 

°bx z (1-v v } D 
x y x 



E M 

_ „ 1 _y 

by (1-v v •. D 
J x y,l y 



a,' = z 



M 



'bxy 



= -zG- 



xy 



The maximum bending and shear stresses occuring at z 
gives 



~2 



M 



X 



a, = ±6 

bX h2 



M 

a... = ±6— ; Tu „„ = +6-X* 



by 



M 
h 2 



"bxy 



( 11 ) 



where and and are given in equations (6) . The 

total stresses are the sums 

a. = a ± a. 
tx x bx 



a. = a ± ct , 
ty y by 



( 12 ) 



T , = T + T, 

txy xy bxy 



2 . 4 Load Boundary Conditions 

The following most general conditions state that the 
edges are subjected to an inplane average compressive load 
per unit length in the y-direction and x-direction of 
magnitudes and respectively. Additionally, all edges 
are subjected to a constant inplane shear load per unit 
length of magnitude S. Referring to Figure 1, 
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at x = ±- 



at y = ± T 



N = S , 
xy 



3 2 F 

3x3y 



= -S 



Inplane load resultant is 



P x = / 



b/2 



3 2 F 



-b/2 3y 



dy = -N b 
2 X 



N = S , 
yx ' 



3 2 F _ - 

3x3y 



Inplane load resultant is 



P = / 

y 



a/2 



3 2 F 



dx = -N a 
-a/2 3x 2 y 



2 . 5 Support Conditions 

Edges simply supported at x = 

Deflection equals zero: w = 0 

External moments equal zero: 



3 2 w . 3 2 w 

+ v 

3x 2 Y 3y 2 



0 



Edges clamped at y = ±j 

Deflection equals zero: w = 0 

rs - _ 

Slope equals zero: = 0 



(13a) 



(13b) 



(14a) 



(14b) 
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3 . Theoretical Analysis 



3 . 1 Analysis Procedure 



The solution method used is identical to that of 
Mansour [9, 10] which is an extension of Levy [4], Coan [3] 
and Yamaki [17]. Briefly, the outline of the procedure is: 



1. Express the deflection of the plate satisfying the 
support boundary conditions, in a double trigono- 
metric series choosing only a finite number of 
terms, w will be a function of unknown nondimen- 
sional coefficients b 

mn 



2. Substitution of this expression into the "compati- 
bility" equation (9) results in a fourth order 
partial-differential equation with the Airy stress 
function F and quadratic functions of the coeffi- 
cient b expressed as coefficients C 

mn ^ pq 



3. A solution for the stress function is assumed which 
satisfies equation (9) and the load boundary con- 
ditions. F will be expressed as functions of 
unknown coefficients d> which in turn are functions 

pq 



of coefficients C 



pq 



To determine the unknown coefficients b , Galerkin's 

mn 

Method is applied to the "equilibrium" equation (7) 



4 . 
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with substitution for w and F by their appropriate 

expressions. Using the orthogonality properties of 

the trigonometric functions, a set of simultaneous, 

non-linear algebraic equations involving cubic 

products of the coefficients b results. 

mn 



3 . 2 Solution 

The deflection of the plate surface can be expressed in 
the form 



Satisfying the boundary conditions (14) , the deflection terms 
take the form 



Differentiation of expression (15) and substitution into 
the "compatibility" equation (9) 



w = h E Z b f (x) g (y) 
m n mn m r n J 



(15) 



V x) 




m = 1 , 3 , 5 




n = 1 , 2 , 3 




xy 3x 2 3y 2 y Sy 4 
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b U £ E E 4b b. .mnij sin— x sin— ? — y 
a 2 b 2 m n 1 j mn ^ a b 



ITT . 2jlT r 

' — x sm — v - 
a b 



sm— -x sm— y - Z I Z Z 4b b. .m 2 j 2 cos^-x 

_ _ • • mn i] J a 

m n l j J 



.,n+l, 2nu , in 2 j tt , 

• I (-1) + cos— g— y] cos— x cos— g— y} 



(16) 



in, i = 1 , 3 , 5 . 

ri/ j = 1/ 2 * 3. 



which can be expressed in the form 



J^ + 2J 
X 3x 4 



3 4 F + j 3 4 F 
y 3y 4 



xy 3x 2 3y 2 



h ir „ „ ^ 2pir 2qiT 

Z Z C cos — — — x cos— y 

a 2 b 2 P q pq a b 



(17) 



p, q = 0 , 1, 2.... 



where C are quadratic functions of the nondimensional 

Ir tl 

coefficients b 

mn 

A particular solution for (17) is assumed in the form 



F = h 2 Z Z <p cos— x cos-^^-y 

p p q r pq a b J 



(18) 



substitution of (18) into equation (17) defines coefficient 



<j) as 

pq 



where 




The stress function F to satisfy the load boundary conditions 
(13) and equation (9) is 



F 



N 



_Z X 2 

2 




Sxy + F p 



( 20 ) 



The method of B.G. Galerkin is applied to the "equilib- 
rium" equation (7) to determine the unknown coefficients b 

mn 

Galerkin' s condition requires that the following equation be 
satisfied by all functions f (x) g (y) 

IT S 



a/2 b/2 

f n f n [D. 



9 4 w 
X 9x 4 



+ 2D 



9 4 w 



xy 9x 2 9y 2 



+ D 



Y 



9 4 w 
9y 4 



- _ ^ IK . 3 2 w 2 9 2 F . 9 2 w _ 9 2 F # 9 2 w 
q ~ 9y 2 * 9x 2 9xSy * 3x9y ” 9x 2 * 9y 2 



( 21 ) 



• f (x) g (y) dxdy = 0 

JL b 



where f (x) g (y) is given be equation (15) . Equation (21) 

3T S 

may also be obtained by applying the principle of virtual 
work. 

Substitution of w and F, expressions (15) and (20) 
respectively, into equation (21) and using the orthogonality 
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properties of the trigonometric functions, there results 



l E [(2s) 
r s 



2 P 4 
0J 

y 



N b 2 

(-*— ) 
tt 2 D 

y 



+ 




N a 2 

(— — ) J b 
* 2d x 



rs 



- ey- & IS b rn {2(-l) n+s r 4 + [r 4 +2p 2 n(2s) 2 r 2 + p 4 (2s) 4 ] 6 ns 



N a 2 
n+s , x 



- 2 »r 2 (-1) “ TS (— — )} + Z I £ £ 4m(i) 2 A (-1) s+1 Ml(.. - S . a . b . .) b 

m r n s tt mr 8J__ _ 2 ~ ran 



ir D 



x 



y tt d. 



xy 



odd 



+ Z Z E 2nm(i) 2 A (-i- + — ) ( 

m r . u mr n+s n-s' 6J 



n+s 



^ H£_ / Sab \ b 

2 r- 1 ran 



y i D 



xy 



+ Z Z E E b {2 (-1) s+1 n 2 r 2 £ <J> 
m n r s mn un+r 



, n 



+ d> + d> ] 

y m-r Y r-m 



~,n 



2 ' 



n 



n+1 



+ 2(-l) il,J 'm 2 s 2 [0, + <J> + <J> ] 

v ' r m+r- r m— r- m 



m+r r m-r 



r r-m 
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( 22 ) 



m, r= 1, 3, 5.... 

n, s 1, 2, 3 . . . . 
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where 6. . is the Kronecker delta, <J> . . = 0 if i < 0 or 

1 J 1 / J 



j < 0, = 0 if n = s and 

J ' n-s 



mr m+r 



if 



m+r 



odd 



or 



mr m-r 



if 



m-r 

2 



odd 



Equation (21) is identical, except for the addition of the 
two shear and two N load terms, to that of Mansour [9] case 
(i) and, if isotropic material properties are substituted, 
reduces to Yamaki [17] case Illb for zero initial curvature. 

Substitution for coefficients $ in equation (22) 
results in a set of simultaneous, non-linear algebraic 
equations involving cubic products of the coefficients b mn . 
Solution of these equations gives the coefficients b and 
hence defines the functions w and F. 



3 . 3 Bending and Membrane Stresses 

The membrane stresses are now determined from equations 
(10) , (18) and (20) . These are 

a = - - h Z £ <p ( “£~“) 2 cos— x cos^3^-y 

x h p q Y pq b a b J 

o = - — h 2 2 4> (^£IL) ^ cos-^— x cos-^^-y 

y h p q T pq a a b J 



xy 



-| - h E £ * (2Ei) (2ai)sin2Elx sin^y (23) 

h p q pq a b a b J 
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The membrane stresses will be represented in the most 
general form, the nondimensional form used by Mansour [9] . 
Using equation (19) 

a 



o* = 



x 



N* 

X 



(^) 

a 2 J 

y 



- Z Z 



12 (1-v v ) 
x y 



C q 2 p 4 

.. P. q . _ 



(24) 



P q 4 (p 4 + 2yp 2 p 2 q 2 + p 4 q 4 ) 



2p7T 2qTT 

■cos — — — X cos — r~ — y 
a b 
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(-^) 



N J 



12(l-v x v y ) 



b 2 J 



x 



and 



(25) 
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YP' 
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p, q = 0, 1, 2.... 



The bending and twisting stresses from equation (11) are 
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(26) 



The normal bending stresses may be expressed as functions of 
nondimensional bending moments M* and M* where 



M 



M* 

x 



X 



D D 



= M* + V 
x_ T 



[_1_2L_X.] 
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— M* 
y y o 
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X o 
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M 
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x 
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M* 
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hD 
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(27) 



(28) 



The bending stresses as functions of the nondimensional bending 



moments are 
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6M* 




D D 

±JLJX. 



6M* hD 




(29) 



The nondimens ional moment due to curvature in the x-direction 

only (M* ) and that due to curvature in the y-direction only 
x o 

(M* ) are determined from (15) and (28) 

Yq 



^ 7T 2 m 2 , mir r , , . n+1 , 2nir , 

M* = E E b cos — x [(-1) + cos-r-y 

x mn 2 nm a 
o p 2 



(30) 



,.j. „ „ . 2 2 i- ITITT 2n7T 

M* = E E 4tt n b cos — x cos— r— y 
y m n mn a b 2 

2 o 



m ■* 1/ 3/ 5 . . . . 
n 1 f 2 , 3 . . . . 



The bending stresses may be expressed as nondimensional bend- 
ing stresses in a form compatible to equations (24) and (25) . 
Using equations (15) and (26) , the nondimensional bending 
and twisting stresses are 



* bx 

at = = + 

bx _ 2 



mir 



I_Ul] 

a 2 J 



± Vm r ^ 2 b cos— x 

2(l-v x v ) m n mn a 



4n 2 v ~ 

[m 2 (-l) n +(m 2 + %-) cos — r~y ] 

e 2 b 



( 31 ) 
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a* = , g . b y . = ± 

by - 2 i 



2 (1-v v ) m n ^mn 
x y 





(32) 




and 




E £ b nm sin 
m n mn 




m = 1 , 3 , 5 



n = 1, 2, 3 



The nondimensional total stresses are therefore, from (12) 



3 . 4 Effective Width 

The effective width of a rectangular plate which has 
buckled is defined as that width of a uniformly stressed 
phantom plate of the same thickness stressed to the same 
maximum stress and sustaining the same total force as the 
real plate. The effective widths a g and b g , in the y and x 





(33) 



T* = T* + T* 

txy xy bxy 
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directions respectively are hence defined as 

N 

a o = a - 77 - 
e ye h 



(34a) 



N 



be = b , 
e xe h 



x 



(34b) 



where a is the edge membrane stress in the y-direction 

(a at x = ± I) and a is the edge membrane stress in the 
y 2 xe ^ 

x-direction (a at y = ± y) . Substitution of (25) evaluated 

X ^ 

£ 

at x = ± y into (34a) gives 
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(35a) 
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substitution of (24) evaluated at y = ± j into (34b) gives 



^e 
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x x y p q 



q 2 (- 1 ) q p 4 C 



pq 
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(35b) 
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or 



N* + 3 (1-v v [E- - -Cq , q + 2 

x x y q 2 

^ q F= 

Clearly, the effective widths a g and are not constant over 
the y and x-directions respectively. Neglecting the small 
amount of change due to the periodic terms in y or x, taking 
only the average values into account [15] , the effective 
widths take the form 



< 3 2 '- 1 »^' C pq , 2p, 

“ — ] COS X 

(p 4 + 2yp z p 2 q 2 + p 4 q 4 ) 
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(36a) 



and 




N* 
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(-l) q 

N* + 3 (1-v v )£ --- - - C „ 

x x y q g2 o , q 



q = 0, 1, 2. 



(36b) 



3 . 5 Principal Stresses 

For the state of plane stress, the maximum and minimum 
principal stresses, and o^r are given by 
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Substituting the nondimensional total stresses, equations 
(33) , the nondimensional principal total stresses are 
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a* 

1,2 



1,2 



b 2 J 

x 



and 



if = 
tx 



tx 



p 4 e 2 



a 



*' = a* 



ty 



ty 



-* » 
txy 



txy 

. yp 2 B 



( 



37 



4. Numerical Solution 



4 . 1 Analytic Solution 

Shultz [15] found that eight deflection terms were 
sufficient to describe a simply supported plate up to an 
aspect ratio of 8. In this solution eight deflection terms 
in the y-direction and one in the x-direction were assumed. 
From equation (15) with m = 1 and n = 1,2,3. ..8, the 
deflection of the plate takes the form 

w = h cos— x [b. . (1+cos^-y) + b, 0 (-1+cos^-y) 
all b 12 b 

+ b 13 (1+cos^-y) + b 14 (-1+cos^y) + b 15 (1+cosij^-y) 

+ b 16 (-1+cos-i^y) + (i +cos Ml y ) + b 18 (-l+cos^~y) 



Similarly, the right hand side of equation (16) is ex- 
panded for m, i = 1; n,j = 1,2,3. ..8. The coefficients C 
are determined by collecting terms on the right hand side of 
equation (16) and matching the coefficients to the series 
on the right hand side of equation (17) for p = 0, 1; 

q=0, 1, 2. ..16. The coefficients X and Y are listed 

q q 

in Table 1 where the coefficients C are given by 

pq 



C = Y + (-l) p X 

pq q q 

p = o, 1 

q=0, 1, 2. ...16 



( 38 ) 
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Substitution of coefficients C into equation (22) 
using equation (19) and with m, r = 1 ; n, s = 1,2,3 .... 8 
produces the following equation after considerable simpli- 
fication 
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( 39 ) 



n,s = 1,2,3 . ...8 
where W = 1.0 if n / s 



W = 2.0 if n = s 
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and o 

x(l - v v ) 

A = 4 x y 

pq (p 4 + 2 yp 2 p 2 q 2 + q 4 p 4 ) 

where 

4 

<p = A C ■£?- 

pq pq pqQJ y 

Equation (39) gives a set of eight simultaneous, non- 
linear equations for the eight coefficients b.^, b^' ^ 13 ' 

b 14' b 15 ' b 16' b 17 ' and b 18 ’ 



4 . 2 Computer Solution 

Equation (39) is programmed on the IBM System 370 model 
165 using Fortran IV level Gl. The subroutine ZEROIN from 
the M.I.T. MATHLIB program library is used to solve the 
system of simultaneous, non-linear equations. This sub- 
program uses an iterative method of solution and convergence 

is achieved when the difference between two successive values 

-12 

of b is less than 0.5x10 
nm 

The foregoing analysis is for the completely general 
case with all the loadings as indicated in figure 1 , however, 
the computer solutions were restricted to certain specific 



cases . 
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Basically, two programs were written to calculate the 
deflection coefficients. The first program produces curves 
of deflection at the center of the plate, effective width 
and bending moment in the y-direction at the center of the 
fixed supports versus inplane edge compressive load N* with 

X 

Q* as a parameter, N* = 0 and S* fixed. Curves of this 
nature, for a range of orthotropic parameters y and ru 
result in a set of design charts. The other program permits 
investigation of the behaviour for any set of loading con- 
ditions with N* = 0. Deflection along the centerline x = 0 
and total principal stresses along the centerlines x = 0 
and y = 0 are plotted. 
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5 . Results 

5 . 1 Design Charts 

The design program was used to generate a set of de- 
sign charts with N* = S* = 0. Virtual plate aspect ratios 
of 1/1.5, 1/2, 1/4 and 1/6 are presented. For each aspect 
ratio a complete set of curves for orthotropic material 
parameters y = g = 1.0 and y = 4.0, g = 0.5 are given. 

Description of design charts: 

5.1.1 Charts of Deflection 

The nondimensional deflection at the center of the 
plate, from equation (15) with x = y = 0, is plotted versus 
the inplane compressive edge load N* with the lateral load 
Q* as a parameter. Figures 2,3,8,9,14,15,20 and 21 give the 
deflection. 

5.1.2 Charts of Effective Width 

The nondimensional effective width given by equation 
(36b) is plotted versus the inplane compressive edge load 
N* with the lateral load Q* as a parameter. Figures 4,5,10, 
11,16,17,22 and 23 give the effective width. 

5.1.3 Charts of Bending Moment 

The nondimensional bending moment M* at the middle of 
the supports y = ± and x = 0, equation (27), is plotted 
versus the inplane compressive edge load N* with the lateral 
load Q* as a parameter. Figures 6,7,12,13,18,19,24 and 25 
give the bending moment. 
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5 . 2 Behaviour Charts 

The behaviour program was used to investigate a few 
general sets of loading conditions with N* = 0 . Virtual 
plate aspect ratios of 1/1.5 and 1/1.25 are presented for 
orthotropic material properties y = g = 1.0. 

Description of behaviour charts: 

5.2.1 Charts of Deflection 

The nondimens ional deflection, from equation (15) , along 

the centerline of the plate (x = 0) is plotted versus y/b 

with either S* or N* as a parameter and Q* fixed. Figures 

X 

26,27,36,37,46, and 47 give the deflection. 

5.2.2 Charts of Maximum and Minimum Total Principal Stresses 
The nondimens ional total principal stresses, equation 

(37) , along the centerlines on the bottom surface of the 
plate (z= -^-) are plotted versus y/b (x = 0) and x/a (y = 0) 
with either S* or N* as a parameter and Q* fixed. Figures 
28-31, 38-41, and 48-51 give the stresses for x = 0 and 
Figures 32-35, 42-45, and 52-55 give the stresses for y = 0. 



c 
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6 . Discussion of Results 

6 . 1 Orthotropic Properties 

Reference [10] contains a complete discussion of the 
rigidity coefficients D^, D^. and and compliance coef- 

ficients J x , J and J X y which are used to express the non- 
dimensional parameters and orthotropic material coeffi- 
cients. Possible approximate formulas for their calculatic 
are given in that reference. The two extreme cases of 
orthotropic plate coefficients are the isotropic case where 
r| = y = 1 . 0 and the grid (intersecting beams, no plate) 

D = 0 and y = 00 • Most plate-stiffener combinations fall 
between n = Y = 1*0 and n = 0.5, y = 4.0; hence the choice 
of parameters for the design charts. 

Comparison of Figures 2 to 25 as to the nondimensionai 
load parameters N* and Q* indicates as p decreases the 
assigned values of N* decrease and Q* increase. N* and Q* 
are indirect functions of the aspect ratio since they are 
nondimensionalized by the plate edges a and b respectively. 
The ranges of values used in the design curves were predi- 
cated on achieving a certain minimum value of effective 
breadth in the calculations. 

6 . 2 Comparison with Existing Solutions 

The first set of design curves to be produced were 
those for the case P=Y = n = l*0> Figures 61,62 and 63, 
for N* with Q* as a parameter and S* = N* = 0. This choice 

X X 

of parameters permitted comparison with the solution obtaine 
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by Mansour [10] . The number of coefficients retained in 
this analysis (only one in the x-direction vice two by 
Mansour) does not favor the solution for the square plate 
case, however, the results of Figures 61, 62 and 63 agree 
to within 3% of the figures 1, 2 and 3 of Mansour [10] . 

The deflection coefficients obtained in reference [10] show 
that the one term assumption only neglects terms that are 
less than 5 1/2% of the primary deflection term b^, which 
accounts for the reasonable agreement in this case. 

For the loading conditions of the design curves, 

Figures 2-23, comparison was made with the solution of 
Levy [5] . Levy calculated deflection at the center of the 
plate and effective width for the isotropic plate of 
3= 1.5, Q*=S*=0 and the same boundary conditions as 
this investigation. Comparison of Levy's results, Figure 3. 
and Table 11 of reference [5] , with the curves for Q* = 0 
from Figures 2 and 3 is graphically displayed in Figures 64 
and 65. Correlation of the two results is very good. 

6 . 3 Design Charts 

6.3.1 Charts of Deflection 

The deflection at the center of the plate is always the 
maximum deflection since only the first deflection mode is 
encountered in these cases. With reference to the curve 
Q* = 0 in the deflection charts, the plate remains unde- 
flected until the critical compressive load N* is reached. 
Increasing N* further buckles the plate and plate deflection 
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increases rapidly. Hence, intersection of the curve Q* = 0 
and the N* axis gives the lowest buckling load. The non- 
linearity of the curves is evident. In the small deflection 
range, ^ a doubling of the lateral load Q* will double 

the deflection, however, in the large deflection range the 
lateral load effect becomes nonlinear with decreasing incre- 
ments of plate deflection for increasing additions of load. 

As p goes from 1/1.5 to 1/6 the magnitude of the maxi- 
mum deflection of the plate is considerably reduced. The 
plate becomes more of a beam supported all around and hence 
"stiffer" than the slightly rectangular plate. The effect 
on the deflection of increasing Q* as p + 1/6 becomes less. 

Comparison of buckling loads and the range of N* values 
over the four aspect ratios shows that for p = 1/1.5 the 
maximum value of N* considered is four times the buckling 
load while it is less than twice the buckling load for 
p = 1/6. The reason for this selection of N* values is 
that the range where the solution is reliable decreases as 
p -* 1/6. Convergence to the correct solution becomes in- 
creasingly more difficult. In addition, the 8 term solution 
begins to become insufficient for p = 1/6 where the eight 
term, b^g, has increased to a value almost 6% of the primary 
term, b^, for the highest set of N* and Q* considered. 

6.3.2 Charts of Effective Width 

By definition, the plate is fully effective in carrying 
the external compressive edge load until it buckles. With 
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reference to the effective width figures, b^/b = 1 for Q* = 

0 until the buckling load is reached. Extensions of the 
curves for Q* / 0 at values of N* less than the buckling 

X 

load are meaningless and should be ignored. After buckling, 
the effectiveness drops off as the plate deflects. Compari- 
son of different aspect ratios shows that the effective 
width decreases more rapidly and takes on significantly 
lower values a p -*• 1/6. 

6.3.3 Charts of Bending Moment 

The bending moment M* at the middle of the fixed sup- 
ports y = ± ^ is calculated because it is expected to be 

large. M* is equal to M* at these points since the curva- 
^ y o 

ture along these supports in the x-direction is zero (M* =0) , 

x o 

hence equation (28) is plotted. With reference to the 
bending moment figures, the bending moment is zero until 
the plate deflects. Increasing Q* increases the bending 
moment in all cases with the nonlinearity associated with 
large deflections. 

6 . 4 Behaviour Charts 

Unfortunately, there are no existing solutions to which 
the results of shear loading may be compared for this par- 
ticular set of boundary conditions. Levy [7,8] and Payer 
[11] considered the simply supported isotropic case. How- 
ever, the results obtained are consistent with those antici- 
pated in the magnitude and shape of the center line 
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deflection curve and values of principal stresses. 

Results are given only for p = 1/1.5 and 1/1.25. 
Investigation of the solution revealed that for virtual 
aspect ratios greater than about 1/1.75, the plate will have 
one symmetric buckle about the diagonal while for smaller* 
values of p antisymmetric buckling shapes are experienced. 

The choice of one term in the x-direction limits the solution 
to some degree in the former case*, however, it cannot begin 
to approximate the true solution in the later case. Hence, 
the narrow range of investigation for p. 

Additionally, with reference to the nondimensionalizing 
of the total stresses, sections 3.3 and 3.5, to present the 
results in a concise and meaningful form, the solution is 
restricted to the isotropic case where y = g = 1.0. This is 
a result of the isotropic aspect ratio 3 appearing expli- 
citly in the nondimensional equations. Nondimensionalizing 
the bending and membrane stresses to eliminate 3 will result 
in a set of incompatable nondimensional stresses which must 
be plotted and evaluated separately, thus increasing the 
number of charts required for each specific case. 

Initially, the total principal stresses at both the top 
and bottom (z = ± j respectively) of the plate were investi- 
gated. Figures 56-60 show the results for p = 0.667, y = n = 
1.0, S* = 14 (above critical load) and N* = Q* = 0. Using 
the Von Mises yield criteria for the biaxial case, o =0, 

*Note: The assumption of the one team solution in the x- 

direction will model a symmetric buckle about the plate 
centerline, not the diagonal. 
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(ag) 2 •= (a*) 2 + (a*) 2 - a* a* (40) 

a combined stress may be calculated to evaluate the relative 
importance of the total stress field at the top or bottom 
of the plate. With reference to Figures 57-60, the table 
of (Uq) Z is as follows: 

y=x=o y =t l- x= ° x=± l- y =0 

Top 12.4 17.7 10.8 

Bottom 12.2 19.5 10.8 



The combined stresses are slightly greater at the top 
than at the bottom for the center of the plate, however, the 
magnitude and difference is considerably greater at the 
fixed supports with the bottom stresses being of higher 
value. By the nature of the upwards (+z) buckling mode, 
this is the expected result. Due to the higher combined 
stress at y = ± j, the behaviour charts were plotted only 
for the principal stresses at the bottom of the plate 




Figures 26-45 show the effect of increasing shear loads 
above the critical load for p = 1/1.5 and 1/1.25. Inplane 
edge loads N* = N* = 0 and cases are presented for both Q* = 

0 and a finite value. Increasing S* and/or Q* increases the 
deflection and total principal stresses throughout the plate. 

Figures 46-55 show the effects of completely general 
sets of loading conditions with increasing inplane edge 
compressive loads N* above the critical load. Comparison 
of the charts for the sub-critical (S* = 6) and super-critical 
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(S* = 14) shear loads clearly indicate the influence of 
shear on the shape of the deflection curve and values of 
principal stress. Increased shear increases deflection and 
principal stresses. In addition., the highest maximum 
principal stresses are observed to occur at the fixed 
support, y = ± -j. 

6 . 5 Examples Demonstrating Use of the Charts 

The following examples are given to demonstrate the use 
of the design and behaviour charts. The approximate formu- 
las for calculation of the rigidity and compliance coeffi- 
cients are discussed in detail in reference [10] . The 
formulas are listed in Table 2. 

6.5.1 Design Example 

Consider the following orthotropic characteristics and 
nondimensional loads for the stiffened plate: 

p = 0.5, y = 4.0, g = 0.5, N * = 4.0, Q*=S*=0 

From Figure 9 

tt 2 D 

The critical load N = 1.8 — 

x 2 

c a 

O 

The center of the plate deflection 

w = 1 . 64h 

where h is the average stiffener's depth plus the 
plate thickness. 

From Figure 11 

The effective width b = 



0.505b 



Edge membrane stress (equations 24,34b and 35b) is 




N* 

x 

12 (1-v 2 ) 



1 4.0 

0.505 *10.9 



a - u . / J 

x e a 2 J 

y 

From Figure 13 

The bending moment at the middle of the fixed support 



M* = -102 

y 



The maximum bending stress at this point 
is 



a 





hD D 

= - 612 -2- 
b 2 hb 2 



(equation 29) 



6.5.2 Behaviour Example 

Consider the following orthotropic characteristics and 
nondimensional loads for an isotropic plate: 

p = 0.667, y = n = 1.0, N* = 3.0, Q* = 8.0, S* = 6.0 

X 

From Figure 46 

The maximum deflection at the center of the plate w = 
1.74 h 

From Figures 48 and 50 

The maximum principal stresses occur at the fixed sup- 
ports y = ± j, x = 0 and are 



H N> 
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and the combined stress 



+ ol - a, a 9 = 54.6 t^- 
212 b 2 J 



X 
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7. Conclusions and Recommendations 



The problem of a large aspect ratio orthotropic plate 
with the short edges fixed and long edges simply supported 
has been modeled for the case of inplane edge compressive 
and lateral loads and subcritical edge shear loads. The 
design charts provide valuable information of use in estab- 
lishing design criteria for the plate-stiffener combination. 

From the behaviour Figures 46 to 55, it is evident that 
the highest principal stresses occur at the fixed supports 
y = ± 2 ’. The design Figures 2 to 25 do not present informa- 
tion to permit the calculation of the membrane stresses at 

the fixed supports. Design charts of a , edge membrane 

^e 

stress in the y-direction (denominator of equation (36a)), 
should be produced. 

As noted in the results, the one term assumption in the 
x-direction will not model the deflection of a plate with 
edge shear as the predominate loading. As the major point 
of this investigation was to model a large aspect ratio 
plate, the solution chosen was an inevitable consequence 
when the large range of p is considered. It is recommended 
that to model a large aspect ratio plate subjected to shear 
loading, close attention must be paid to the important 
deflection terms that should be included in the solution. 
Levy [8] confined himself to a M 1.5 and used 14 selected 
terms in the deflection equation. 
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Additionally, the eight term solution is just suffi- 
cient to model the case for the orthotropic plate aspect 
ratio of 1/6. 
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TABLE 1 - COEFFICIENTS C 

pq 




V q + (-1)% 



p = 0,1 

q = 0 , 1 , 2 ... 16 



where X and Y are given by 

q q 

X„ = b 2 + 4b? 0 + 9b 2 + 16b 2 + 25b 2 + 36b 2 + 49b 2 + 64b 2 

o 11 12 13 14 15 16 17 18 

Y = -X 
o o 



X 1 4b ll b 12 + 12b 12 b 13 + 24b 13 b 14 + 40b 14 b 15 + 60b 15 b 16 
+ 84b 16 b 1? + H2b 17 b 18 

Y 1 = 2b ll ( “ b ll ~ 2 ° 12 ~ b 13 + b 14 " b 15 + b 16 “ b 17 + b 18 ) 

“ 13b 12 b 13 " 25b 13 b 14 “ 41b 14 b 15 " 61b 15 b l6 ” 85b 16 b 17 “ 113b 17 b 18 



X^ = 



-b ll + 6b ll b 13 + 16b 12 b 14 + 30b 13 b 15 + 48b 14 b 16 



+ 70b 15 b 17 + 96b 16 b 18 



Y 2 8b 12 ( " b ll + b 12 b 13 2 b 14 b 15 + b 16 " b 17 + b 18 ) 



” b l! 10b ll b 13 34b 13 b 15 52b 14 b 16 74b 15 b 17 



lOObigbig 



X 3 _4b ll b 12 + 8b ll b 14 + 20b 12 b 15 + 36b 13 b 16 + 56b 14 b 17 



+ 80b 1 5 b 18 
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Y 3 18b 13 ( “ b ll + b 12 “ b 13 +b 14 “ b 15 “ 2 b 16 b 17 +b 18 ) 

_5b ll b 12 “ 17b ll b 14 “ 29b 12 b 15 ~ 65b 14 b 17 “ 89b 15 b 18 



X, = 



6b ll b 13 + 10b ll b 15 ” 4b 12 + 24b 12 b 16 + 42b 13 b 17 



+ 64b 14 b 18 



Y, = 



3 '* b 14 (- b n + b 12 b 13 + b 14 - b 15 + b 16 - b 1? - 2 b i8 ) 



-10b b , - 26b 11 b lc - 4b* - 40b 1o b 1(; . - 58b-._b.._ 

11 13 11 15 12 12 16 13 17 



x 

Ul 

II 


+ 
i — i 

■Q 

i—i 
i— i 

A 

00 

i 


12b ll b 16 


12b 12 b 13 


+ 28b 12 b 1? + 48b 13 b 18 


Y 5 = 


5 °bi 5 (-bn 


+ b 12 “ b 13 


1 

i — 1 

A 

+ 


b 15 + b 16 " b 17 + b 18 } 




• 17b ll b 14 - 


37b ll b 16 


13b 12 b 13 


53b 12 b 17 - 73b i3 b i8 



X 6 -10b ll b 15 + 14b ll b 17 “ 16b 12 b 14 + 32b 12 b 18 ~ 9b 13 
Y 6 = 72b 16 ( “ b ll + b 12 “ b 13 + b 14 “ b 15 + b 16 " b 17 + b 18* 



-26b i;L b 15 - 50b 11 b 17 - 20b 12 b 14 68b i2 b i8 “ 9b 13 



X 7 = " 12b ll b 16 + 16b ll b 18 ‘ 20b 12 b 15 “ 24b 13 b 14 
Y 7 = 98b 17 (-b ll + b 12 “ b 13 + b 14 “ b 15 + b 16 “ b 17 + b 18* 
- 37b 11 b 16 - 65b 11 b 18 “ 29b 12 b 15 “ 25b 13 b 14 



x 8 -14b n b i 7 “ 24b l2 b 16 “ 30b 13 b 15 ” 16b 14 
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Y 8 128b 18 ( “ b ll + b 12 “ b 13 + b 14 “ b 15 + b 16 “ b 17 + b 18* 
- 5 ° b H b 17 - 40b 12 b 16 - 34b 13 b 15 - 16b^ 



X 9 


=-16 b n b ig - 


28b 12 b 17 


■ 36b 13 b 16 • 


40b 14 b 15 


Y 9 


= - 65b ll b 18 ' 


53b 12 b 17 


’ 45b 13 b 16 ‘ 


41b 14 b 15 


O 

i — 1 

X 


= ' 32b 12 b 18 ‘ 


42b 13 b 17 


- 48b 14 b 16 - 


25b 15 


Y 10 


=_68b 12 b 18 ‘ 


58b 13 b 17 


- 52b l4 b 16 ‘ 


25b 15 


rH 

i — 1 

X 


=-48bi3 b i8 - 


56b 14 b 17 


- 60 b 15 b 16 




1 — 1 
rH 
>i 


= ~ 73b 13 b 18 


65b 14 b 17 


- 61t, 15 b 16 




X 12 


= ’ 64b 14 b 18 ‘ 


70b 15 b 17 


- 36b 16 




Y 12 


=-80 b i4 b lg 


74b 15 b 17 


' 36b 16 




X 13 


= - 80b 15 b 18 - 


84b 16 b 17 






CO 

rH 

>H 


= - 89b 15 b 18 ' 


85b 16 b 17 






i — 1 

X 


= - 96b 16 b 18 " 


49b 17 






Y 14 


=- 100b 16 b 18 ■ 


- 49b i7 







X 15 “ 112b 17 b 18 

Y 15 = “ 113b 17 b 18 
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x is = - 64b « 
X 16 = X 16 
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TABLE 2 - RIGIDITY AND COMPLIANCE COEFFICIENTS 



D - 
x 



El 

x 

S (1-v 2 ) 

X 



El 



D 



y S (1-v 2 ) 

y 



p “ C- 



I S 

JL 
I s 

x y 



1 



D 



n = 



xy 



.1 



D D ' 
x y 



I I 

pxpy 



I I 
x y 



J - 



x Eh 



J - 



y Eh 



x 



xy 



1 r l+v 
E 1 h 






P h 



Y = 



xy 



- (1 + v) ^ 



h K 
x y 



- v- 



h h 

x y 



J J 

x y 



2h h 
h = 



h + h 
x y 



61 



TABLE 3 - INTEGRALS 



b/2 a/2 

f / 

0 0 



r r miT ru , , ab 

/ / cos — x cos — x dxdy = -5- 

a a 8 



= 0 



m = r 
m ^ r 



b/2 a/2 

f / 

0 0 



r , mir r tt 2 stt , , 

/ / cos — x cos — x cos— r— y dxdy 

« a a b 



= 0 



b/2 a/2 



/ / cos— x cos— x cos ^v "y dxdy 

0 0 a a b 



rir 



2niT 



= 0 



b/2 a/2 

r / 

0 0 



r r rnTT it tt 2nir 2stt , , ab 

/ / cos— x cos— x cos — g — y cos -g - y dxdy = > m=r, n=s 



= 0 m^r or n^s 



/ / cos— x cos— x cos— x cos ^- 11 y dxdy = 0 

0 0 a a a b * 



b/2 a/2 

/ / 

0 0 



rTT 2sn , - _ 

cos— x cos y dxdy = 0 



b/2 a/2 

/ / 

0 0 



r-1 

ru , , ab § -i \ 2 
cos— X dxdy = 27F<-1> 



= 0 



r odd 
otherwise 



, ^ ^ • mTT 2piT r it 2qTT 2 stt , 

/ / cos — x cos — =- — x cos — x cos— y cos— r— y dxdy 

„ a a a b J b- 2 ^ 



0 0 



ab , „ 

= * m+r = 2p, q = s 

= ^>|m-r| = 2p , q=s, m > r 



ab 

32 



>|r-m| = 2 p, q=s, r > m 



0 otherwise 
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b/2 a/2 m7r 2 p 7r 

/ f cos — x cos — — — x 

a a 

0 0 



rir 2mr ^air , , 

cos— x cos-—— y cos-g— y dxdy 

ab 



2qTr_ 



32 

ab 

32 

ab 

32 

= 0 



m+r=2p , q = n 
m - r | = 2p, q = n, m > r 

> [r - m| = 2p, q = n, r > m 
otherwise 



b/2 a/2 . _ 

/ / sin^-x sin— — x cos— x sin — P— - y sin— 3^y dxdy 

oo a a a b " 



ab 

32 



m+r=2p , n=q 



= ' | m - r| = 2p , n=q, m > r 

= |r - m| = 2p, n=q, r > m 



= 0 



otherwise 



e t mu 2 ptt r-rr 2nir 2 q 7 j 2 stj , , 

/ / cos — x cos-- — x cos— x cos — i— — -y cos - -g— y cos- ? -~ y dxdy 

n n a a a b J b ^ b J J 



0 0 



_ v 

= -gj > m + r = 2p , n + s = q 

ab , „ ■ - 

= -rj > m + r = 2p, n - s I = q, n > s 



64 

ab 

64 

ab 



j m + r = 2p, I s - n|= q, s > n 



- •£-£> | m - r|= 2p, n + s = q, m>r 



ab 

64' 


II 

u 

1 

g 


2p t 


|n - 


s l 


II 


m>r, n>s 


ab 
64 ' 


1 m - r | = 


2p , 


|s - 


n| 


II 


m>r, s>n 


ab 

" 64' 


1 

_3_ 

II 


2p f 


n + 


s 


= q/ 


r > m 
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= .jr^ j |r - to | = 2p, |n - s|= q, r>m, n>s 



= 

64 
= 0 



/ |r - m|= 2p,|s - n|= q, r>m, s>n 
otherwise 



b/2 a/2 > , 2pu 

J J O ~i r\ c? l r\ £. 

0 0 



,m . 2mi . 2qn 



2stt. 



sin— x sm— x cos— x sin— g— y sin-^— y cos— g— y dxdy 

ab 

> m+r=2p,|n-s|=q,n>s 



) m + r = 2p, n + s = q 

ab 
64 

ab 



64 

ab 

64 

ab 

64 

ab 

“64 

ab 

“64 

ab 

“64 

ab 

64 

= 0 



) m + r = 2p,|s -n|=q, s > n 



m-r 



= 2p, n+s=q,m>r 



m - r | = 


2p, 


l n 


- 


s l 


1 = 


q. 


m>r , 


n>s 


3 

i 

h 

II 


2p t 


I s 


- 


n l 




q/ 


m>r , 


s>n 


h 

I 

_3_ 

II 


2p , 


n 


+ 


s 


= 


q. 


r > 


m 


1 

_3_ 

11 


2p , 


In 


- 


s l 




q / 


r>m, 


n> s 


h 

1 

_3_ 

11 


2p, 


I s 


- 


n l 




q/ 


r>m. 


s>n 



otherwise 



b/2 a/2 

/ / 

0 0 



. imr rir . 2nu , , 

sm — x cos — x sm— r— y dxdy 
a a b 2 J 

= ab f 1 ^ 

2 W } ' 



nn 

ab 



<=l>, 



mr 
= 0 



2 m-r 



m+r 

2 

m-r 



odd , n odd 
odd, n odd 



otherwise 
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b/2 a//. ' m7r r 7r _ 2n7T 2 sir , , 

J J sm — x cos — -x sin-T— y cos-r-y dxdy 

0 0 a a D D 



ab 1 . , 1 1 , m+r 

0 2 l m+r 'n+s n-s^ J 2 

2 7T 



odd, n+s odd 



. ab , 1 . , 1 + 1 . 

2tt 2 ' n+s n-s^ > 



m-r 



odd, n+s odd 



= 0 otherwise 
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Table 4 • 


- Computer Symbols 


Fortran Symbol 


Variable 


B(N) 


b. 

l,n 


RHO 


P 


GAMMA 


Y 


ETA 


n 


QSTAR 


Q* 


NXSTAR 


N* 

X 


NYSTAR 


N* 

Y 


SSTAR 


S* 


W/H 


w/h 


AE/A 


a e /a 


BE/B 


b /b 
e' 


MYSTAR 


M* 

y 


SIGMA STAR 


a* or a* 


CPQ 


c pq 


APQ 


A pq 


RANG 


Max. value of N* or S* 



RANGQ 



Max. value of Q* 
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Figure 1 



p-\ s: tri s: 



3.2-, 
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N* 

x 



Figures 2 & 3 
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2-0 4.5 7.0 S . 5 12.0 




N* 

x 



Figures 4 & 5 



M* 

y 







N x 



Figures 6 & 7 






3.2 
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0 2.0 4.0 6.0 8.0 




N* 

x 



Figures 8 & 9 
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N* 

x 



Figures 10 & 11 



20 

-60 

140 

220 

300 

10 

-70 

150 

230 

310 



p = 0.50 




2.0 4.0 6.0 8.0 



N* 

x 



Figures 12 & 13 



31 S *** 
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N* 

x 



Figures 14 & 15 



1.02 
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b 

e 

b 



i 





N* 

x 






Figures 16 & 17 



20 
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Figures 18 & 19 



zr\z, 
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N x 



Figures 20 & 21 
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K 



Figures 22 & 23 



20 

160 

340 

520 

700 

300 

50 

200 

450 

700 
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P 

Y 

n 

Q* 



0.167 

1.0 

1.0 

0,50,100,150 



.0 



T 



0.5 1.0 



1.5 2.0 




P 

Y 

ri 

Q* 



0.167 

4.0 

1.0 

0,50,100,150 



0.5 



1.0 



1.5 



2.0 



N* 

x 



Figures 24 & 25 



1 .00 
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p = 0.667 y = p = 1.0 N* = 0 




Figures 26 & 27 



SfOH SIRR STGMR STRR 
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p = 0.667 y = n = 1.0 N* = 0 





Figures 28 & 29 



St.GMri STRR 5IGMR STRR 

5.0 7.5 10.C .0 7.5 3 .0 
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p = 0.667 • y = n = 1.0 N* = 0 

o 





Figures 30 & 31 



STGMR STAR STGMA STAR 

-2.5 -1.7 -.8 .1 r 3.2 -2.4 -1.6 
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p = 0.667 y = n = 1.0 N* = 0 





Figures 32 & 33 



yuis tiwois ybis bwois 
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p = 0.667 y = n = 1.0 N* = 0 

rvj 




Figures 34 & 35 



1 .05 1 .40 
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p = 0.80 y = n = 1.0 N* = 0 




Figures 36 & 37 



SIGMA STAR SIGMA STAR 

-2.2 -L.H -.6 .2 - 2.1 -2.0 -1.3 
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p = 0.80 y = r) = i.o N* = 0 





Figures 38 & 39 



SIGMfl STAR SIGMA STAR 
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p = 0.80 y = 0 “ 1.0 N* = 0 




Figures 40 & 41 



SIGMA STAR SIGMA STAR 
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p = 0.80 Y = n = 1 . 0 N*=0 





Figures 42 & 43 



SIGMA STAR SIGMA STAR 
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p = 0.80 Y = n = 1.0 N* = 0 




Figures 44 & 45 



H/M M/M 



p = 0.667 



Y = n = 1.0 
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N* 

y 



0 




Figures 46 & 47 



SIGMA STAR SIGMA STAR 

-5.0 -3.7 -2.4 -l.l .2 r 3.U -2.2 -1.4 
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p = 0.667 y = n = 1.0 N* = 0 





Figures 48 & 49 



SIGMA STAR STGMA STAR 

6. a 9.5 13.0 -l .0 ] .5 11 .Cl 
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p = 0.667 y = n = 1.0 N*=0 

y 





Figures 50 & 51 



SIGMA STAR SIGMA STAR 
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p = 0.667 



y = n = l.o 



N* = 0 

y 





Figures 52 & 53 



SIGMA STAR STGMA STAR 
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p = 0.667 y = n = 1.0 N* = 0 

C3 
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p = 0.667 



Y = n = 1. 0 Q*=N*=0 

S* = 14 




Figure 56 



J 



SIGMH STRR SIGMA STAR 

2 .0 3 .0 4 .0 -4 .u -3 -U -2 .1) 
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p = 0.667 y = n = 1.0 Q* = N* = 0 




- .33 



.15 



~ 



r/B 



T 
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I 

.52 



Figures 57 & 58 



[Gm STRR SIGMR ST^R 
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p = 0.667 y = n = 1.0 q* = n* = 0 




Figures 59 & 60 



3*13 
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M* 

y 
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-20 



-40 



-60 



-80 




N* 



Figures 61 & 62 




1.02 



0.89 



0.76 



0.63 



0.50 




N* 

y 



Figure 63 
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* x 




SU: 



■i 



Figure 64 
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d) 



A 



I* 



Figure 65 
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9 . Appendices 



APPENDIX A 



Details of Derivation of Basic Equations of 
Large Deflections 
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Definition of External Forces and Moments 
Shearing force per unit length parallel to x axis 



Q x - 


h/2 

/ T dz 

-h/2 XZ 



Shearing force per unit length parallel to y axis 



o 

II 


h/2 

/ T dz 

-h/2 yZ 



Bending moments per unit length of sections of a plate per 
pendicular to x and y axis respectively 



II 


h/2 

/ z a dz 
-h/2 x 


II 


h/2 

/ za dz 

-h/2 y 



Twisting moment per unit length of section of a plate per- 
pendicular to x and y axis respectively 



M 

xy 


h/2 

= -/ zx dz 

-h/2 Xy 


M 

yx 


h/2 

= -/ ZT dz 

-h/2 yX 



i 



104 



Strain - Displacement Relations 



Consideration of strain in the middle of the plate 




For large deflections, the effect of deflection on the 
strains has to be included in the strain-displacement re- 
lations. The accompanying figure shows a typical elongation 
of a linear element of length dx due to displacement in the 
x and z directions from bending. 

Elongation = dl - ^dx + i(|^) 2 dx 

o X 4 oX 



The strain in the x direction is therefore 



e 



x 



3u . 1 /3W\ 2 
9x 2 v 9x ; 



(a) 



The strain in the y direction is likewise 



- ~ f 2 

y 3y 2 v 3y ; 



(b) 



The shearing strain due to displacements u and v is 



3u 9 v 
3y 9x 



It can be shown that the shearing strain due to displacement 



w is 
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9w a 8v/ 
3x * 8y 



The total shear strain in the middle plane is 

v * 2* + 9X + # gw 

T xy gy gx gx gy 

For small angles the displacements u and v can be represented 
by the change in displacement w as 



Substitution into (a) , (b) and (c) gives 



x 




(2)* 



e 



y 





*Numbered equations correspond to those used in the main 
text. 



Equilibrium of Forces and Moments 




The above diagrams are superimposed to make one loading 
condition. For equilibrium in the z-direction: 
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EF = 0 
z 



3Q 

x 



3Q. 



g-^dxdy + ■ dydx + q dxdy = 0 



!^x 

3x 



3Q 

= _ 

9y 



(a) 



Moment equilibrium about the x-axis 



EM =0 
xx 



3M. 



xy 



3M 



3Q. 



3x 



dxdy - -r-^dydx + Q dxdy + -r-^dydxdy + qdxdyO (dy) 
o y y d y 



= 0 



Neglecting the last two terms since they are of higher order 
gives 

3M 3M 

- - - x ^ - -X- + Q = 0 (b) 

3 x 3 y y 

Moment equilibrium about the y-axis and similarly neglecting 
small terms 



EM. 



yy 



o 

3y ^ 3x w x 



0 

0 



2. Forces in the middle plane of the plate (Membrane Stresses) 
The forces acting in the middle of the plate may have a 
considerable effect on the bending of the plate and must be 
considered. 
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Assuming no body forces or tangential forces 

ZF = 0 
x 

3N 3N 

-^dxdy + -^dydx = 0 



ZF = 0 

y 



3N 3N 

— - ^ -dydx + ■ ■■ x ^ dxdy = 0 
ay a x 



and from symmetry noting N = N 

J r xy yx 



i 



i 

i 
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9N 
x 

9x 



+ 



9N 




9N 



9Y 



9N 

. *y 

9x 



(3) 



Equations (3) are independent of equations (a) , (b) and (c) 

and can be treated separately. 

Projecting normal forces onto the z axis and taking 
into account the bending of the plate and resulting small 
angles 



-N dy|^ + (N + . 

x J 9 x x 9 x 



x , w gw , 9 2 w , . , 
dx) (f- + - 2 — dx) dy 

SX 9x 2 



Similarly projecting normal forces N onto the z axis 

9N 

-N dx-— + (N + — - ^dy) (— - + - - - - dy) dx 



9y 



9y 



9y 



9y- 



Neglecting higher than second order terms gives 

^ 2 9N 

M 9 w, , , x 9w, , 

N — -dxdy + -» — • 7 T— dxdy 
Xg x 2 9x 9x 2 

~ 2 9N ~ 

N — ^dydx + • -^dydx 

y 9y2 y 9y 9y y 



(d) 



(e) 



Treating the projections of the shearing forces onto the 
z axis the same way 



( 1 ) for N with slope of deflection surface in the y 



direction on the two opposite sides of the element as 

dx 



9w , 9w , 9 2 w 
TTy an 9y 9x9y 
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N 



9 2 w 



9N 



xy9x9y 



dxdy + iF 



3 w , 



(2) for N with slope of deflection surface in the x 

direction on the two opposite sides of the element as 

<ay 

9N 



3w , 9w , 9 2 w 
3x an 3x 3x3y 



N 



9 2 w 



yx3x9y 



dydx + 



yx 



3Y 



9w, , 



combining (1) and (2) gives 

9N 



2N 



9 2 w 
xy9x9y 



dxdy + 



& + 



9N 



xy 9w, 
9y 9x 



xdy 



(f) 



Differentiation of (b) with respect to y and of (c) with 
respect to x and substitution into (a) eliminates shearing 



forces Q and Q 
x y 



9 2 M 



x 



9 2 M. 



yx 



9 2 M 



y _ 



9 2 M. 



xy 



9x‘ 



3x9y 



9y‘ 



3x3y 



q 



Noting that M = - M , adding expressions (d) , (e) and (f) 

yx xy 

to the load qdxdy originally defined and making use of (3) 



9 2 M, 



x 



3 2 M 3 2 M 

2 + — ac 



9x‘ 



9x3y 



9y‘ 



= -(q + N 



9 2 w 

x^ 2 
9x z 



+ 2N 



9 2 w 
xy9x9y 



+ N 

y 9y 2 



(4) 
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Moment - Curvature Relations 



Substitution of equations (2) for and y in equations 

Cl) gives* 



o 



x 



x 



1-v V 

X y 



I-z 



9 2 w 
9x 2 



+ 



1 /3w> 2 

2 l 9x ; 



zv 



9 2 w 
r 3y 2 



1 ,9W. 2 i 

2 V y<37> ] 



0 = 
Y 



1-v V 

X y 



[-z 



9 2 



w 



9 Y‘ 



+ IA 2 - 

2 l 9y ; 



zv 



9 2 



w 



X 9x 2 



1 ,9W. 2i 

2 v x ( sr> 1 



(a) 



T 

xy 



-2zG 



9 2 w 
9x9y 



+ G 



9w 9w 
9x 9y 



Taking the moments as defined before and integrating 



M 



x 



h/2 

/ Z0 dz 
-h/2 x 




h/2 

/ zo dz 
-h/2 Y 




h/2 

-/ ZT 
-h/2 




Even functions will drop out over the integration interval 
leaving the same expressions for moments M^, M and as 
for a plate undergoing pure bending and only small deflec- 
tions. The result is 

*Numbers in parenthesis correspond to equations in the main 
text. 
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M = 
x 



M = 

y 



- D x ( 



•D ( 

y 



9 2 w + 


3 2 w. 




v , 


9x 2 


y 3y 2 


9 2 w 


9 2 w, 


4- 


v 


3y 2 


9x 2 



( 6 ) 



M = 2C^_ 
xy 9 xgy 



where the rigicity coefficients and D y are defined as 



E h 3 
x 



D x - 12(l-v V ) 
x y 



D = 



E y h 3 



y " 12 (1-v v ) 
J x y 



and 



C = 



Gh 3 

12 



"Equilibrium 11 Equation 

Substitution of equations (6) into (4) gives 



D + 2D 



9**w + D 3 ** w 



x_,i, xy, v2 ,„ 2 y„„4 



9x 



9x^9y' 



9y" 



= q + N 



ll™. + 2N + N_ 3 2 W 



k 9x‘ 



xygxgy y 3y2 



2D = D v + D v +4C 
xy x y y x 



where 
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Substitution of the stress function defined by equations (5) 
^ 3 4 w 



3 4 w 



X 3x** 



+ 2D — + D 

xy 3x 2 3y 2 y 



9 4 w 
3y 4 



_ - , 3 2 F 3 2 w 0 3 2 F 
= q + - 2 

3y 2 3x 2 



3 2 



w 



3x3y 3x3y 



3 Z F 3 2 w 

— - - * - 

3x 2 3y 2 



(7) 



"Compatibility” Equation 



The forces N , N , and N in the middle plane of the 
x y xy ^ 

plate depend on the strain due to bending as well as the 
external forces applied in the xy plane. 

Again assuming no body forces and requiring load q is 
perpendicular, equations (3) apply for equilibrium in the 
middle xy plane 



3N 



x 

3x 



+ 



3N 



xy 



3y 



0 





(3) 



The corresponding strain components are those of equations (2) 

_ _ + 1 2 
x 3x 2 v 9x / 



e 



y 




1 ,3w. 2 
2 ^3y ; 



( 2 ) 



_ 9u + 3v^ 3w # 3w 

Y xy 3y 3x 3x * 3y 



Differentiating equations (2) and combining to eliminate 
u and v results in 
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3 2 e 



3y ' 



8 2 Y 



xy 



3 2 e 



3x3y 



3x 



_y = ( 9 w ) 2 _ 

2 v 3x3y' 



3 2 w 3 2 w 



3x ' 



9y : 



( 8 ) 



Using Hooke's Law to relate strain and forces N 



N 



N 



e = 



- v 



'x hE y hE 

x J y 



e = 



N 

He 



N 



- V 



X 



x hE 



x 



(a) 



'xy 



N 

= x y 

hG 



Differentiation of equations (a) and substitution into 



equation (8) with use of equation (3) 



3 2 N 



x 3x 2 



3 2 N 



+ 2 J 



xy 



xy 3x3y 



+ J 

y 



3 2 N 
x 

3y 2 



= ( 



3 2 w } 
3x3y ; 



3 2 w 3 2 w 

■ — - • ■ — - 

3x 2 3y 2 



where 





1 

E h 
x 



and 



2 J 



xy 



_1_ 

Gh 



- v J - 

x y 



v J 
y x 



Substitution of the stress function defined by equations (5) 

ji!Z + 2 j a-*- +j rg. (!!«)» (9) 

x 3x' xy 3x 2 8y 2 y 3y" 3x3y 3x 2 3y 2 
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Effective Width 





Figure b 



Consider the simple case of a rectangular plate simply 
supported on all edges where the loaded edge remains straight 
in the plane of the plate at all times. If the load is below 
the buckling load, the stresses will be distributed evenly 
as shown in figure a. For post-buckling loads the center 
of the plate will exhibit less compressive strain than the 
edges because of large deflections of the center of the 
plate. The stress distribution is shown in figure b. The 
effective width relates the maximum stress o uniformly 
distributed along a phantom plate sustaining the same total 
load as the real plate, thicknesses being the same. 

Therefore 



a a 
e max 



a a 



Where a is the average edge stress. For this particular case 
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N 



-?r- = 5 and 



a max °ye 



Where is the edge membrane stress in the y-direction 

(a at x = ± ^) • Substitution into the above relation gives 
the result 



N 

y 

a a = a ~ 
e ye h 



(34a) 



Similarly, for loads in the x-direction 

N 

bo = b — 
e xe h 



(34b) 
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1. Determination of Coefficients C 
pq 

Substitute deflection expression (37) into the right 
hand side of equation (16) . The right hand side of equation 
(16) becomes* 



Term #1 




. • 8n . . IOt: . . 12tt 

+ 4b,sin-r-y + 5b c sin-T— -y + 6b,sin-r— y 
4 b J 5 b J 6 b J 




Term #2 






+4b2Cos-^-y + Sb^cos— y + 16b^cos~y + 25b^cos-~^y 



+36b 6 cos^^y + 49b^cos-^^-y + 64bgCos^jj2y] } 



(a) 



*Note: since m = 1 for all cases, for convenience b 

represented by b . 



l,n 



is 



n 
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Each of the two terms of equation (a) is multiplied out and 
like terms are collected in the form cos^p-y for 
q = 0,1, 2... 16 using the function product relations 

4sinasin$ = 2 cos(a~3) - 2cos(a+$) 

4cosacos8 = 2 cos(a-|3) + 2cos(a+6) 



The first and second terms produce coefficients 2X and 2Y 

q q 

respectively. Using the function relationships 



. 2 IT 1 1 2 7T 

Sin 4 — X = yr - vCOS X 

3. Z Z 3 



cos 2 ^x = i + icos^-x 

3 Z Z 3 



there results 

First term: 



X cos— 3— y 
q b 



„ 2q7T 2 7T 

X cos— r— y cos — x 
q b J a 



Second term: 



Y cos^32-y + Y cos— 3— y cos— x 



Comparison with the right hand side of equation (17) , the 
coefficients C are 

pq 



C = Y + X 
o,q q q 



C. = Y - X 

i,q q q 



are 
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or 




Y g + (-l)?X q 



P = 0/1 
q = 0,1,2. 



16 



2 . Application of the Method off BY G. Galerkin to Equation 
(7) 

Method of B.G. Galerkin (from reference [18]) 

It is required to determine the solution of the equation 

L(w) = 0 

where L is a differential operator in two variables, which 
solution satisfies homogeneous boundary conditions. 

The approximate solution sought is in the form 

i 

w (x,y) = I b n F n (x,y) 
n=l 

where F n (x,y) (n = l,2...i) is a system of functions chosen 
to satisfy the boundary conditions and b^ are undetermined 
coefficients. Consider the functions F n (x,y) to be linearly 
independent over a given region (no one of the functions can 
be expressed as a linear combination of the others) . For 
w(x,y) to be the exact solution of the given equation, it is 
necessary that 



L(w) = 0 



122 



If L (w) is continuous, this requires that the expression 
L (w) is orthogonal to all functions of the system F n (x,y) 

(n = l,2...i). With i constants, b^, t^**** 3 !' i conditions 
of orthogonality can be satisfied. 

Requiring that the two functions are orthogonal over a 
region produces the following system of equations 



// L(w)F (x,y) dxdy 
R 

= ff L( £ b n F n (x,y) )F n (x,y)dxdy = 0 
R n 

(n = 1 , 2 . . . i) 

from which the coefficients b n can be determined. These 
coefficients will define the solution w(x,y) . 

The method of B.G. Galerkin can also be obtained from 
the principle of virtual work. 

Application of the method of Galerkin to equation (7) 
requires that the following equation be satisfied by all 
functions f (x)g (y) 

2T 2 



a/2 

/ 

0 



b/2 

/ 

0 




3 **w 
3x 4 



+ 2D 



xy 



3 4 w 
3x 2 3y 



2 



+ 




q 



3 2 F 

3y 2 




+ 2 



3 2 F 
3x3y 



3 2 w 
3x3y 



3 2 F 3 2 W j 
3x 2 3y 2 



X f (x) g (y) dxdy = 

L 2d 



0 



( 21 ) 
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Substitution of w and F, expressions (15) and (20) 
respectively 



a/2 b/2 

t 0 ?§ (hD x =os^y]cosHLx 



, <^* v rt ri / 2mr \ 2 /Ttitt » 2 v 2 n tt 

+ 2 ZZ ) z (— ) z b_ cos— r— y cos — x 



xy mn b 



mn 



_nnr_ 

a 



v vv ri ri / 2 n TT v 4 1 2 TT TT ITT TT 

+ h D Z Z (-r— ) 4 b cos — r- — y cos— x - q 
y m n b mn b a 



v 3 v* v* v v /HITT \ 2 / 2 q TT \ 2 , t r /• n \ n^ - 1 . 2nTT ^ 

- h Z Z Z Z ( — ) (--—) <b b [ (-1) + cos-r-Y 

m n p q a b y pq mn 7 b 1 J 



.2qTT 



mir 



.2PIL 



• cos -' y cos — x cos—* — x 
b J a a 



- , v v ,2mT. 2i_ 2mr mir 

N h E Z (■- , ) z b cos— r— y cos — x 

ymn b mn b a 



h 3 E E (^^-) 2 (I P— ) 2 (j) b cos^^-y cos — if— -y cos^x cos-^-x 



m n b a y pq mn 



. 2qiT. 

b 



. irni . 

a 



■IPJL, 

a 



Ct u v v /ITlTT \ 2 r / i \ nH“ 1 . 2nTT i ITlTT 

N h Z Z ( — ) z b [(-1) + cos — r" — y ] cos — x 

xmn a mn b J a 



o o v* _ , • 2nTT % ,imr \ i • 2mT • ititt 

- 2Sh S, 5 <—> ( -i-> b mn sln — S' sln T* 

t 2 h' [ [ [ I ( 2 Bi) (— ) ( 221 ) (231) 4 , b sin 2 lly sin 2 gly 
mnpq b a a b r pq mn b 2 b 2 



. ITITT 

sin — x sm 

a a 



2pTT *, r, ,vS+L 2STT t rTT . , n 

— ^x} • [ (-1) + cos— r— -yjcos — x dxdy = 0 

a Da 



m, r = 1,3,5. .. . 

n , s 1,2,3.... 
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The orthogonality properties of the trigonometric 
functions must be used to evaluate the integral of equation 
(21) . By employing the function relations 

1 i 

cosA + cosB = 2cos-j(A + B)cos-j(A - B) 
sinA + sinB = 2sin-j(A + B)cos^(A - B) 



The trigonometric products may be expressed in typical forms 
a/2 

/ cos— x cos— x dx 

0 a a 

l a/2 

= -rf I2cos— (m+r)x + cos— (m-r)x + cos— (r-m) x] dx 
4 q a a a 



and 



a/2 

, . mir rTT n 

f sm — x cos — x dx 

0 a a 



, a/2 

= 4 / [2sin— (m+r)x + sin— (m-r)x 
4' q a a 



sin— (r-m) x] dx 
a 



where integration is carried out for all m and r. By 
combining terms that integrate to similar values, these 
integrals are equivalent to 
a/2 

[2cos^-|m+r | x + 2cos^-|m-r|x + 2cos^-| r-m| x] dx 

and 



, a/2 

jj-f I2sin^-|m+r | x + 2sin^-|m-r|x - 2sin^-| r-m|x]dx 



where the second two terms in each expression integrate to 
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zero if r > m or m > r respectively. Manipulation of the 
trigonometric products on the right hand side of equation 
(21) into these typical relations permits the reduction of 
the equation into the convenient form of equation (22) . 
Table 3 lists all the specific integrals from equation (21) 
with their respective values. 
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COMPUTER PROGRAMS 
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The computer program to solve the simultaneous equations 
(39) is written in Fortran IV Level G1 for the IBM system 
370. The program consists of the main and three subprograms — 
MATHLIB subroutine ZEROIN, and subroutines EVAL and CPQ. 

Table 4 lists the important computer symbols used. 

ZEROIN is a Fortran IV subprogram, from the M.I.T. 
mathematical library, which computes the solution vector 
b of a set of N simultaneous, nonlinear equations 
N 

EF n (b) = 0 using double-precision arithmetic. b is obtained 

by an iterative method beginning with estimated values of 
the solution vector and iteration is performed as 

b k+1 = b k + D k 

1c k 1c lc 

with the vector D = J *F(b ) where F (b ) represents the 

k. k 

vector of function values at the point b and J represents 
an approximation to the Jacobian at b . Convergence is 
tested with the expression 

| b k+1 - b k | 2 S 0.5 x 10“ 1 2 | b k+1 | 2 

Subroutine EVAL generates the eight simultaneous 
equations and calculates the vector of function values for 
a given solution vector b. EVAL is called repeatedly by 



subroutine ZEROIN. 
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Subroutine CPQ is in turn called by subroutine EVAL, 
with a solution vector b, to calculate the coefficients C 

pq 

that are used in the simultaneous equations. 

For the starting values of the solution vector b, the 
main program uses a well guessed set of values set in the 
program. If the lateral load is not zero (Q* ^ 0) (see 
note 1) or the present loading conditions did not change 
significantly from the last set, the previously computed 
solution vector b is used for the starting values. 

If ZEROIN cannot achieve convergence in 50 iterations, 
the main program alters the initial set of starting values 
and returns them to ZEROIN. Normally convergence is obtained 
within the first 50 iterations, however, never more that a 
second set of values is required. 



Part 1 - Design Charts 

The main program varies the nondimensional loads Q* 
and N* or S* for a particular set of plate parameters. 

Tables of nondimensional deflection coefficients, nondimen- 
sional deflection at the center of the plate, nondimensional 
effective width and nondimensional bending moment (y- 
direction) at the middle of the supports (y = ± b/2) are 
outputs. Additionally, CALCOMP plots of the last three 
tables are produced. 
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Part 2 - Behaviour Programs 

The main program reads the input parameters for 4 sets 
of loading conditions on a particular plate and outputs 
CALCOMP plots of nondimensional plate deflection along the 
centerline (x = 0) and nondimensional total principal 
stresses along the centerlines (x = 0 and y = 0) . These 
programs are completely general and can be used with any 
set of loading conditions. 

The design programs must be used with care. With 
reference to equations (36a) and (36b) , the effective width 
is only defined for N* ^ 0 or N* ^ 0 respectively. Hence, 
variation of S* with N* = 0 will give a value of zero for 
the effective width. The design programs are written so 
that no computation or plot of effective width will be 
made for any variation of S* with N* fixed. Additionally, 
deflection at the center of the plate and the y-direction 
bending moment may not be significant parameters in this 
case. It is recommended that the "behaviour" program be 
used if variation of the shear load, S*, is desired. 
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Program Descriptions 



Main Program 1 - Design (N* = 0) 

Input one card per set of complete design curves set up as 
follows : 

Column Parameter 



1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

*Note: Either N* or S* may 

same run. 



RHO 

GAMMA 

ETA 

*NXSTAR (Enter Negative Value 
to vary) 

*SSTAR (Enter negative value 
to vary) 

RANG 

RANGQ 

be varied , but not both in the 



Output 

1. Tables 

a. b-i ^ 

1 1 n 

b. ^ at y = x = 0 

c. b (only if N* is varied) 

0 X 

d. M* at y = ± -j, x = 0 

y ^ 



i 
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2. CALCOMP Plots 

Tables b, c and d above. 

3. Values 

N* or S* = 0 to RANG in 10 increments of RANG/20 and 
Q* = 0, RANGQ/6, RANGQ/3 , RANGQ 

Termination 

A blank card as the last data card is required. 
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Main Program 2 - Design (N* = 0) 

Input one card per set of complete design curves set up as 
follows : 

Column Parameter 



1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

*Note: Either N* or S* may 

same run. 



RHO 

GAMMA 

ETA 

*NYSTAR (Enter negative value 
to vary) 

*SSTAR (Enter negative value 
to vary) 

RANG 

RANGQ 

be varied, but not both in the 



Output 

1. Tables 

a. b. 

1 r n 

2. CALCOMP Plots 

a. ^ at y = x = 0 

b. a (only if N* is varied) 

e y 

c. M* at y = i x = 0 
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3. Values 

N* or S* = 0 to RANG in 10 increments of RANG/20 and 
Q* = 0, RANGQ/6 , RANGQ/3 , RANGQ 

Termination 

A blank card as the last data card is required. 
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Main Program 3 - Behaviour 

Input four cards per program run, each set up as follows: 



Column 




Parameter 


1-10 


RHO 


) 


11-20 


GAMMA 


Imust be same 


21-30 


ETA j 


1 


31-40 


NSTAR 


(see note 2) 


41-50 


SSTAR 




51-60 


QSTAR 




61-70 


*TB 





*Note: If TB other than zero, stresses at top of plate will 

be computed. TB = zero defaults to bottom of plate. 

Output 

1. Tables 

a. t>. 

1 ,n 

b. | at y = x = 0 

c. cj and a* at y = x = 0 

2. CALCOMP Plots 

w n - 1 . y .,1 

a. K at x = 0 for - ^ S § < + j 

b. c* and ot at x = 0 for - i 5 ^ 1 + i- 

1 Z A cl A 
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c. a* and cr£ at y = 0 for 



1 < * < + 1 

2 ~ a 2 



Termination 

A blank card as the fifth card is required. 
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Subroutine EVAL 

Two subroutine programs EVAL - one for N* = 0 , the 

other for N* = 0. The appropriate EVAL program must be 
used with the corresponding main programs. 



Subroutine ZEROIN & CPQ 

Completely general subprograms. Use with all combina- 
tions of main programs and EVAL. 



NOTE 1 

If the lateral load is equal to zero (Q* = 0) a possible 

but trivial solution is that all the coefficients b =0. 

mn 

Care must be exercised to avoid using b =0 for starting 
values in this case. Additionally, under certain circum- 
stances the solution readily converges to zero rather than 
the appropriate non-trivial solution. The programmer should 
be aware of this possibility. 



NOTE 2 

The behaviour program may be used with either N* or N* 

y x 

as a loading, but not both. The main program must be used 
with the appropriate subroutine EVAL. 



DES IGN PROGRAM 
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MAIN PROGRAM 1 







- DESIGN 
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NO PLOT IF A SOLUTION DIO NOT CONVERGE 
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20 CONTI NUE 
N0= 1 1-MM 

ARRANGE DATA TO BE PLOTTED INTO VECTORS 
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13 CONTINUE 

WRITE (KO, 250 ) 

INITIALIZE LOADING CONSTANTS 
QSTAP = 0.D0 
VARY CSTAR 



VARY=0.n0 

IF (NST.LT. 0.0) NSTAR=VARY 
IF(SST.LT.O.O) SSTAR=VARY 
VARY NSTAR OR SSTAR 
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T ( 3) = 0 . DO 
T ( 4) = J.D3 
T ( 7 ) =0 .no 
F S=F ( S ) 
S2=2.D0*FS**2 



ST2=( 2 .DO*FS ) **2 
SL = (-1.P0)*«( S + l) 

CALL CPO(S, D, CIS, COS) 

TC42 P = (APQ(G1 t FS ) * CIS +2 .DO^-APQ (GO, FS)*COS )*S2 
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